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Abstract 

m ; 

^ . We study the scaling regime of nucleon - anti-nucleon annihilation 

x/^ \ into a deeply virtual photon and a meson, pN — >■ 7*7r, in the forward 

\ kinematics, where \t\ « ~ s. We obtain the leading twist ampli- 

^ I tude in the kinematical region where it factorizes into an antiproton 

■ distribution amplitude, a short-distance matrix element related to 

^ . nucleon form factor and the long-distance dominated transition dis- 
tribution amplitudes which describe the nucleon to meson transition. 

, We give the evolution equation for these transition distribution 

Q \ amplitudes. The scaling of the cross section of this process may be 

>-Ch ■ tested at the proposed GSI intense anti-proton beam facility FAIR 

• with the PANDA or PAX detectors. We comment on related pro- 

^ . cesses such as vrA^ — > A'7* and 7*A^ — > N'n which may be exper- 

^ I imentally studied at intense meson beams facilities and at JLab or 

■ ■ ■ ' Hermes, respectively. 



1 Introduction 

In a recent paper [IJ, we have shown that factorization theorems ^ for 
exclusive processes apply to the case of the reaction 7r+ — > 7* 7 in the 
kinematical regime where the virtual photon is highly virtual (of the order 
of the energy squared of the reaction) but the momentum transfer t is small. 
We also advocated the extension of this approach to the reaction pp 7*7 
and to virtual Compton scattering in backwards kinematics. This enlarges 
the successful description of deep exclusive reactions in terms of distribution 
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Figure 1 : The factorization of the annihilation process p p ^ -f* tt into the an- 
tiproton distribution amplitude (DA), the hard subprocess amplitude (Th) and a 
baryon — > meson transition distribution amplitude (TDA) . 



amplitudes (DA) [3] and/or generalized parton distributions (GPD) [H [5] 
on the one side and perturbatively calculable coefficient functions describing 
hard scattering at the partonic level on the other side. We want here to 
describe along the same lines the reaction 

Pik)N{p)-^^*{q)n{p') (1) 

in the near forward region and for large virtual photon invariant mass Q, 
which may be studied in detail at GSI [6|. Such an extension of the GPD 
framework has already been advocated in the pioneering work of [7j. 



In Ref. [T] , we defined the tt — >^ 7 leading twist transition distribution ampli- 
tudes (TDAs) from the matrix elements 



=0,z=0 



(2) 



where the Wilson line [y; z] = F exp ig(y — z) J^dt n^A^{ty + (1 — t)z) 

provides the QCD-gauge invariance for non local operators and equals unity 
in a light-like (axial) gauge. In a similar way, we shall define in section 2 the 
nucleon to meson TDAs from the matrix elements 



(vrl g"(zi) [zi] Zq] q'^{z2) \z2] Zq] q^{z^) [z^] Zq] \p) 



=0, zj =0 



(3) 



The pN 7*7r amplitude at small momentum transfer is then propor- 
tional to the TDAs T{xi,^,t), where (i=l,2,3) denote the light cone mo- 
mentum fractions carried by participant quarks, and ^ is the skewedness 
parameter connected with xb hy 

2, — Xb 

in the Bjorken hmit. It reads schematically 

M{Q\ ^,t)^ J dxdy(t>(yi, Q'')T„{xi, y,, g')r(x,, t, Q^) , (5) 

where (j){yi, Q^) is the antiproton distribution amplitude and Th the hard 
scattering amplitude, calculated in the colinear approximation. We shall 
show in section 3 that these TDAs obey QCD evolution equations, which, 
as always, follow from the renormalization group equation of an appropriate 
operator, in our case of the three quark operator. Their dependence is thus 
completely under control. We calculate in section 4 the hard amplitude and 
derive some phenomenological model-independent predictions of our picture. 
In section 5, we comment on processes related by crossing, such as ttN — > 
and 7*A^ — ^> N'n which may be experimentally studied at intense meson 
beams facilities and at JLab or Hermes, respectively. 



2 The N ^ TV transition distribution ampli- 
tude 

Let us take a closer look at the transition distribution amplitudes from a 
nucleon to a pseudoscalar meson. A similar description of the antiproton to 
meson TDA may be straightforwardly deduced from our study. For their defi- 
nition we introduce light-cone coordinates v"^ — ± v^) / \pl and transverse 
components Vt = for any four-vector v. The skewedness variable 

^ = — with A = p' and P = {p + p')/2 describes the loss of plus- 

momentum of the incident hadron in the proton — > meson transition. We 
parametrize the quark momenta as shown on Fig. 1. The fractions of -|- mo- 
menta are labelled Xi, X2 and X3, and their supports are within [— 
Momentum conservation implies (we restrict to the case ^ > ) : 

5^x, = 2e. (6) 
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The fields with positive momentum fractions, > 0, describe creation of 
quarks, whereas those with negative momentum fractions, Xj < 0, the ab- 
sorption of antiquarks. The spinorial and Lorentz decomposition of the ma- 
trix element follows the same line as in the case of the baryon distribution 
amplitude [9l |10]. Because of that let us first recall the definition of the 
proton DA at leading twist 



4(0|e*^-^<(zi nyp{z2 n)d^^{z, n)\B{p, s)) (7) 

where z, j, k are color indices and n is the light cone + direction. The vector 
is a light-like vectors {rt^ = 0) which together with p^^ defines the light-cone 
kinematics. 

We then define the leading twist TDAs for the p — * 7r° transition as : 



4{7r'{p')\e'='ul{z,n)uj,{z,n)d'^{z,n) \p{p, s)) (8) 

^ z+=0,ZT=0 

In 



2/. 
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M 



/■y ; 



where a'^'^ = i/2['~f^,'~f'^], C is the charge conjugation matrix and B the nu- 
cleon spinor. P = P^7^, the vector A = p' — p has - in the massless limit - 
the transverse components 

= (g^"' - -^(P^n" + P"n^))A, . 
Pn 

/tt is the pion decay constant ( = 93 MeV) and is the constant which 
determines the value of the nucleon wave function at the origin, and which 
has been estimated through QCD sum rules to be of order 5.3 ■ 10~^ GeV^ 
[TT| . The first three terms in ([8]) are the only ones surviving the forward limit 
At — > 0. The constants in front of these three terms have been chosen in 
reference to the soft pion limit results (see below). With these conventions 
each function V{ziP ■ n), A{ziP ■ n), T{ziP ■ n) is dimensionless. Finally 
let us note that the number of leading twist TDAs in ([H]) corresponds to 
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eight independent helicity amplitudes related to the matrix element in 
(272 = 8). 

For the n ^ TDA the analogous expression has the form 



2+=0, 2t=0 



4(vr {p')\e'^''u'^{zin)u^p{z2n)d''^{z^n) \n{p,s)) 
Jn_ 

+T3-(PV^,C)„;3(a^^A^fi)^ + ^(A^PV^,C),;3(AtP)-, . (9) 

One might reexpress and Vj~ (respectively and A", respectively T° 
and T~) in terms of the isospin 1/2 and 3/2 quantities V^~^^'^ and V^^"^''^ (re- 
spectively A\~^^'^ and y4[~^^^, respectively t/"^''^ and t/"^''^). Simple isospin 
rotation enables to deduce the TDAs for p tt^ and n ^ 7r° transitions from 
Eqs. (El) and ([9]). We do not write down the p 7r~ transition (which is 
pure isospin 3/2 exchange ) since it does not contribute to the process under 
study. 

Each TDA can then be Fourier transformed to get the usual representa- 
tion in terms of the momentum fractions, through the relation 



F{ziP ■ n) = j £x5{xi + X2 + X3- 2^)e-^^'^^^'"> F{xi, ^ (10) 

-1+5 

where F stands for Vi, Ai, Ti and / d'^x = J dxidx2dx^5{2S, — xi — X2 — 0:3). 



3 Evolution equations 

QCD radiative corrections lead as usual to logarithmic scaling violations. 
The scale dependence of the proton to meson TDAs is governed by evolution 
equations which are an extension of the evolution equations for usual DAs 
and GPDs [3],|3l|5]. The derivation of evolution equation for TDAs proceeds 
in an analogous way as for DAs and for GPDs therefore we sketch only 
essential steps. 
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The non-local three-quark operators relevant for TDAs and their evolu- 
tion, see [TU] for details and notation which we follow, involve quark fields 
having definite chirality or helicity 

gT(^) = i(l±7^)g. (11) 

The separation of "minus" components of quark fields [9] leading to the 
dominant twist-2 contribution is achieved by the substitution q —>■ hq, with 
n = n^j^. There are two relevant operators in our problem: the first one 
corresponds to the case where the three quarks have total helicity 1/2 

B'S^izi, z,, zs) = e^^''=(ng|)„(2;in)(ngj)^(z2n)(ngl,)7(^3ri) (12) 
and the second one with the total helicity 3/2 

Ba^'yi^i, Z2, zs) = e'^'' (nqj) ^{ziu) (nqj) p{z2n) {nql)^{z3n) . (13) 

For simplicity, we have assumed that all quarks in f[T^ . f[T51) have different 
flavours. The conditions imposed by flavour symmetry do not influence the 
evolution equations but lead to certain symmetry properties of TDAs. Since 
operators (fT2l) . (fT3|l belong to different representations of the Lorentz group 
they do not mix with each other. 

The operators B, ( JT2l) and ( JT3ll . satisfy the renormalisation group equa- 
tion 

l^-^B = H-B (14) 
a jj, 

with H being an integral [10] operatoi0 

cv 

H = --^ [(1 + l/iVj n + 3C^/2] (15) 
Ztt 

in which the second term ~ corresponds to the self-energy corrections of 
each quark fleld. The operator Ti. acts in different way on the operator (fT3!) 
and on the operator (fT2!) . In the flrst case it is determined by contributions 
from one loop Feynman diagrams describing in the Feynman gauge the vertex 
corrections corresponding to gluon exchanges between quark flelds and gluons 
forming Wilson lines as 

^3/2 = ^12 + ^23 + ^13' (16) 
^We restored in Eq. US]) the factor -Q;s/(27r) absent in Eq. (2.23) of [10]. 
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where 



1 

n\^B{zi) = - I — {a[B{z%Z2,z^)-B{z^,Z2,Z:,)] 

J OL 


+a [5(^1, ^2^1, ^3) -5(^1,^2,^3)]} , (17) 

with a = 1 — a, z^^ = ZiO. + ZkO.. In the case of ( fT2l) the operator Ti is 
determined not only by above contributions but also by those ones which 
correspond to Feynman diagrams with gluon exchange between quark lines 
having opposite chiralities (i.e. in our case between lines (1,2) and (2,3)) and 
it can be written as 



where 



with 



7^1/2=7^3/2-7^12-^2 3' (18) 

ni^B{z,) = / i?(zn,^2"i,^3) , (19) 

1 

j Va = j (iai (ia2 'i«3 (5(1 — — tt2 — ^3) • (20) 


From the renormalisation group equation (fT4l) and (fTS!) we derive the 
corresponding equation for the matrix element of operators B between states 
relevant for the process under study. The subsequent use of parametrisations 
of these matrix elements in terms of different TDAs, as those given by Eqs. ([8]) 
or ([9]), results in the evolution equations for TDAs Vi, Ai, Ti. 

As a definite example we present the evolution equation for the set of 
TDAs denoted as F^^\xi) which are related to the matrix element of three 
quark operator with total helicity 1/2 

(7rO(p')|e^^''(ngJ)a(^in)(ngj)^(;^2n)(ngl),(^3n)|Ar(p,s)) 

being parametrised according to Eq. ([H]). The evolution equation has the 
form 
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-1+5 

1+5 



xip{x[,xi) 
x[{x'i — Xi) 



1+5 



-1+5 



X2P{X2,X2) 

4(4 - 




' Xip{x[,Xi) 


+ j dx'^ 
^ -1+5 


'X3P(X'3, 




x[{x'i — Xi) 




3^3). 


X2P{X2,X2) 
_x'2ix2 -X2)_ 


1+5 


"X3p(x'3, 


3^3) " 
3^3). 



F^^\x^^,X2,x'^ 



FT-'-'^(xi,X2,X3) 



-1+5 



1+5 



1+5 ^ -1+5 
1+5 1+5 

dX2^p(x2,X2) + / ciX3^p(x3,X3) I FT-'^^(xi,X2,X3) 

-1+5 -1+5 



1 



, 3^3 



The integration region in each integral is restricted in two ways. Firstly, the 
support of integrands is defined by functions p(x, y) = 9{x > y > 0) — 9{x < 
y < 0), with 6{x > y > 0) = 6{x > y)0{y > 0). This function p{x,y) 
is a generalization of analogous one which appears in equations describing 
the pure evolution ERBL [4J. The second condition is the requirement that 
although not denoted as a variable of integration the variables x- must satisfy 
the condition that x'^ G [— 1 + ^, 1 + e.g. in the first integral over x[ on 
the rhs of fl21l) the variable = 2^ must belong to the interval 

x^G [-1+^,1 

The evolution equation for the set of TDAs F^^^{xi), which correspond 
to the case where three quarks have total helicity 3/2, is obtained from (j2T|) 
by neglecting two last lines. 

The results of this evolution are different in the various Xi sectors. In 
particular, when all Xi > one is in the same kinematics as the usual ERBL 
equation for the baryons, with the simple Xi — > Xi/2^ rescaling. The solutions 
of the Eq. (I2T!) in this ERBL region are thus well known and are expressed 
in terms of Appell polynomials P„(xj/2^): 



F(Xi,^,p^) = XiX2X35(xi+X2+X3-2^) ^0„P„(Xi/2^) 



as(Po) 



In /bo 



(22) 
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where 60 = H/SAi'c — 2/3?2/, 7„ are the corresponding anomalous dimensions 
and 0„(/io) are dimensionless nonperturbative parameters. 

When one or two Xi < , the solutions of evolution equation fl2T]) are 
unknown: they will deserve further study. Although asymptotic solutions are 
simple, one should not use them without caution in phenomenological studies, 
since it is known that the asymptotic solution of the proton distribution 
amplitude (namely, the DA proportional to X1X2X3 like the first term in fl221) ) 
does not allow for a good description of form factors at accessible values of 
Q^. Thus one should not insist on the asymptotic form but instead use 
nonperturbative techniques such as QCD sum rules or lattice calculations to 
get boundary values to insert in the evolution equations, and maybe solve 
them by applying some methods based on conformal symmetry used in [TU] 
or by methods proposed in [12j. 



4 Hard amplitude and cross section estimates 

Since the leading twist antiproton distribution amplitude selects the helicity 
±1 state for the hard scattered quarks, and since the photon coupling does 
not modify the helicity of these quarks, the three quarks extracted from the 
proton by the TDA have also a total helicity of ±|. Moreover, in this first 
study, we concentrate on terms which are not vanishing at zero A^, i.e. 
to the contributions of the three first TDAs in Eq. ([HD, namely Vi,Ai,Ti. 
The hard amplitude is then at leading twist straightforwardly deduced from 
the studies of proton form factors at high Q^. At leading order in a^, the 
amplitude Ai^ for the reaction 

p{k,X)p{p, s) l*{q)TT^{p') 

may be read off from baryonic form factor calculations ^ il5j as 

= -tepFiQ\^,t)v{k,\)YlMp,s) 

F{Q\^.t) = f^ ^^"%fP^^ j AI:T.(x.,,,) (23) 



1+5 



with 



4 v + r 



3(2e-Xi+«e)2(x3 + «e)(l-l/i)2i/3 
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T2 



4 r 

3(xi +ie)(2^-X2 + ie)(x3 + ie)?/i(l - ?/2)z/3 

4 V 

3 (xi + ie){2^ -X3 + ie){x3 + ie)?/i(l - yi)y3 
_4 V 

3 {x2 + ie){2^ -X3 + ie){x3 + ie)y2{l - yi)y3 



^2 V /I 1 



3 (2^ -X3 + iey{l - 7/3)2 \ (xi + ie)?/i (x2 + ie)y2 

T = 2 V + T 

6 3(2^_3;^+,e)2(x2 + ^e)(l-i/i)2y2 ^ ^ 



T7 

Tg 
Tg 
Tio 

with 



2 v + r 



3 (2^ - xi + ie)2 


(x2 + ie)(l - ?/i)2?/2 




1 


V 




3 (xi + ie)(x2 + 


2e)(2e-X3 + «e)i/i(l- 


1/1)2/2 


1 


r 




3(xi +ie)(2e- 


- xi + ze)(x2 + «e)yi(l 


- ^2)2/2 


1 


V 





3 (xi + 2e)(x2 + ie){2^ - xi + ie)x3t/it/2(l - 2/3) 



V(x,,i/„e,t) = [Viy,^-A{yi)] ■ [V^ixj,^,t) - A,ix„^,t)] , (25) 

r(x„y,,e,t) = -12[T(|/,)].[Ti(x,-,e,t)] , (26) 

and / d^y = J dyidy2dy35{l — 2/1 — 2/2 — 2/3)- Note that the electromagnetic 
vector current has been replaced in the amplitude (l23l) by an axial vector 
current, due to the presence of the outgoing pseudo scalar meson vr-meson in 
the TDA dHD. 

A model independent result of our analysis is the scaling law for the 
amplitude : 

valid up to logarithmic corrections due to DA and TDA anomalous dimen- 
sions. On the other hand the ratio : 



M{Q\0-^^. (27) 



da{pp ^ l+l-7r°)/dQ'^ 
da{pp in-)/dQ'^ ' 
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should be almost independent. 

Another interesting consequence of our framework is the dominance of 
the transverse polarization of the virtual photon, which results in a specific 
angular distribution of the lepton pair in its rest frame, namely 

adu 

where 6 is the usual emission angle of the lepton in the virtual photon center 
of mass frame. 

Since the electromagnetic form factor Ff is obtained with the same ex- 
pressions as Eqs. (!23|) . with the replacement of the axial vector current 
by vector one and with 2/^ 1, and T — > 4[T(?/j)] • [T(a;j)] one may obtain 
an estimate of the threshold cross section in terms of the electromagnetic 
form factor , once a reasonable form of V, A and T are chosen. We shall not 
enter more the phenomenology in this letter but will study it in a forthcoming 
work. 



5 Related processes and conclusions 

We have defined the new transition distribution amplitudes N —>■ n, i.e. 
which parametrize the matrix elements of light-cone operators between a 
baryon and a meson states; this generalizes the concept of GPDs for non- 
diagonal transitions. Similar matrix element was introduced also in Ref. [7] 
in discussion of the exclusive production of forward baryons off nucleons. 
Obviously related processes are 

• Firstly, the exclusive lepton pair production 

T^N iV'7* , (29) 

in the kinematical regime where the outgoing nucleon is almost colinear 
to the incoming meson, which is the backward region of the reaction 
ttN — > 7*A^' studied in Ref. [13], and which may be studied in an intense 
pion beam facility such as the project JPark. 

• Secondly, the same framework may be applied to 

'y*N N'tt , (30) 
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in the kinematical regime where the outgoing meson is almost cohnear 
to the incoming proton, which is the spacehke analog of the previous 
case. 

• Thirdly, the crossed [t — > s) version of the N 7c TDA describes the 
exclusive fragmentation of three colinear quarks in a Baryon-meson 
pair. Analogously to the two meson generalized distribution amplitude 
[1], it is an interesting tool to access reactions where an emerging iso- 
lated baryon is replaced by a baryon-meson system. An exemple of its 
usefulness is the calculation of the impact factor of a hard diffractive 
reaction (at large values of t) where a baryon projectile is transformed 
into a baryon-meson ejectile. In this way, the crossed TDA describes 
the partonic content of a continuum or resonating baryon-meson state. 

Many other processes may be studied in the same way, where the tt meson 
is replaced with other mesons, or /and the outgoing nucleon is replaced with 
other baryons. 

The introduction of transition distribution amplitudes thus allows to 
study the reactions 

pp — > l^l^n^ pn l^l^n^ 

along the same lines as the timelike electromagnetic baryon form factors. 
The applicability of a perturbative QCD approach to the proton form fac- 
tor at accessible energies has been controversial for years |T3]. We pretend 
here neither that experimental data have shown that the perturbative ap- 
proach is successful, nor that next to leading order corrections will succeed 
in describing them. We just want to point out that a new phenomenology 
of related but different reactions is now possible to try to understand the 
present puzzle. Moreover, the study of form factors has been the subject of 
many developments after the pioneering papers [Sj [15]. In particular, the 
calculation of next to leading logarithm corrections [16], the proposal of an 
optimal renormalization scale fixing [T7], the importance of soft gluon re- 
summation [IH] and the discussion of the timelike vs spacelike aspects [T^ . 
have put to a firmer basis the QCD description of these objects(for a short 
review, see [2Dj). These interesting developments ought to be applied now 
to the meson production processes. This necessary effort should help a phe- 
nomenological analysis to be more constrained so that one may clearly see if 
a (quite sophisticated) perturbative QCD approach is indeed relevant to the 
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kinematical range which may be accessible at GSI-FAIR, with ^/s ~ 32GeV^ 
in the target rest frame and ^/s ^ 207GeV^ in the accelerator mode. On the 
other hand, an alternative more phenomenological point of view has been 
also developped to describe form factors in terms of generalized parton dis- 
tributions and the so-called " Feynman mechanism" |21] . This framework has 
recently been enlarged [22] to describe reactions ([T]) at fixed angle through 
generalized distribution amplitudes [23] . 
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